Temperature dependence of the collective mode and its influence on the band 

splitting in bilayer cuprates 
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The recently observed bilayer splitting in high-Tc cuprates is analyzed within a model where the 
charge carriers are coupled to a phenomenological bosonic spectrum which interpolates between the 
marginal Fermi liquid structure and collective mode type behavior as a function of temperature. 
We argue that the origin of the collective mode is associated with dynamic incommensurate charge 
density waves. Moreover it is shown that the resulting temperature dependence of the self-energy 
E is in good agreement with that extracted from angle-resolved photoemission data. 
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I. INTRODUCTION 

One of the long standing unsolved problems in angle- 
resolved photoemission spectroscopy (ARPES) of high- 
Tc cuprates was the detection of band doubling in bi- 
layer materials. In these compounds one naturally ex- 
pects a splitting of the quasiparticle (QP) dispersion in 
an antibonding (AB) and bonding band (BB) due to fi- 
nite hopping of charge carriers between the two Cu02 
planes. Only recently this feature has been measured in 
BiaSraCaCuaOg+s (Bi2212)ii2i^ above and below T^. The 
splitting turns out to be essentially zero along the diag- 
onals (0, 0) — > (tt, tt) and (in the normal state) aquires a 
value between 88meV^ and llQmeV^ at the (7r,0) point 
of the Brillouin zone (BZ). In addition it seems that the 
magnitude of the energy splitting is constant over a large 
range of doping."* Most interestingly it turns out that 
for the same doping the splitting is only 22meV below 
T(<i contrary to the expectation that the hopping between 
planes tj_ should be essentially temperature independent. 

In this paper we show that the apparent reduction 
of the bilayer splitting below Tc can be explained by 
the assumption that the charge carriers are coupled to 
bosons which frequency distribution evolves from a fea- 
tureless structure well above Tc (i.e. approaching that 
of a marginal Fermi liquid^) to a sharp, mode-like ex- 
citation at small temperatures. As a consequence, low 
energy QP excitations are essentially confined to ener- 
gies up to the mode frequency ujq which thus constitutes 
an upper limit for the splitting between AB and BB QP 
peaks (see Fig. 4). Moreover, the damping of excitations 
for energies larger than loq leads to a broad hump which 
turns out to be mainly determined by the bonding band 
(cf. Fig. 3c). This scenario is thus in agreement with 
recent ARPES experimental where this peak-dip hump 
structure at (tt, 0) has been resolved below Tc in addition 
to the bilayer splitting. 

Numerous experiments now support the existence of 
a collective mode in the high-Tc materials most promi- 
nently exemplified by the appearance of a kink in the elec- 
tron dispersion along (0, 0) (tt, tt) as seen in ARPESSii 
(cf. also Fig. 2). Concerning the origin of the mode the 



various proposals include a magnetic resonance (see e.g. 
Ref. 0), the coupling to phonons^ or incommensurate 
charge-density waves (ICDW)^ as the source of the asso- 
ciated scattering. It has been argued^ that the difference 
between peak and dip position of the photoemission line 
shape at (tt, 0) decreases as a function of doping as it 
is expected for the magnetic resonance frequency. An 
analagous analysis^ came to the conclusion that the fre- 
quency is essentially constant (if not slightly decreasing) 
upon doping thus resembling the behavior of a phonon 
mode. We have recently fittedii ARPES data taking into 
account also the feedback of the doping dependent su- 
perconducting (SC) gap. The extracted mode frequency 
is a strongly decreasing function of hole concentration 
reaching a value of w ^ 40meV. . .50meV around optimal 
doping. This behavior is thus consistent with a quantum 
critical point (QCP) located near optimal doping where 
the associated order parameter corresponds to ICDW 
formation.**^ In fact, the formation of stripe order is now 
well established in a large class of cuprate compounds 
(for an overview see e.g. Ref. IT3I) . Most interestingly the 
coexistence of charge modulations and superconductivity 
has recently also been reported for the Bi2212 compound 
using scanning tunneling microscopy. ^^'^'"^ These experi- 
ments strongly support our view that ICDW scattering 
can also be responsible for the mode-like features as seen 
in ARPES experiments. In the following we substanti- 
ate this scenario by a detailed analysis of the temper- 
ature dependent collective mode which allows us to fit 
the experimentally observed bilayer splitting in the su- 
perconducting and normal state respectively by simply 
tuning the broadening of our phenomoneological bosonic 
spectrum. 



The paper is organized as follows. In Sec. II we ana- 
lyze available experimental data in order to elucidate the 
frequency and temperature dependence of the collective 
mode in Bi2212 materials. Based on these results the 
temperature dependent bilayer splitting is investigated 
within the ICDW mode scenario in Sec. III. We finally 
conclude our discussion in Sec. IV. 
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II. TEMPERATURE DEPENDENCE OF THE 
COLLECTIVE MODE 

The coupling of charge carriers to a collective mode 
is reflected in a drop of the imaginary part of the elec- 
tronic self-energy Im E below some characteristic mode 
frequency luq- In Refs. ITMitI the ARPES hue width, 
which is directly related to Im S, has been analyzed for 
an optimally doped Bi2212 sample (Tc=90K) as a func- 
tion of frequency and temperature. In fact, it was ob- 
served that in the SC state Im E is strongly reduced be- 
low 70...80meV which compares rather well with the drop 
of the scattering rate as obtained from infrared reflectiv- 
ity measurements,-'^^ Moreover, with increasing temper- 
ature there is a gradual change of Im E towards a lin- 
ear frequency dependence above Tc. The corresponding 
marginal Fermi liquid (MFL) structure of the ARPES 
line shape is also discussed in Ref. 

In order to interpolate between these two limiting self- 
energies we consider the coupling of charge carriers to the 
following bosonic spectrum: 



(1) 



where JLj(,_r(w) is a normalized distribution function cen- 
tered around uiq with halfwidth P. Assuming particle- 
hole symmetry with respect to the chemical potential 
and neglecting for the moment the influence of the SC 
gap^, evaluation of the T — self-energy yields 



ImE(w) = X^N{0) 



dyF^^^r{v) (2) 



ReE(w) = -A^iV(O) / dvF^^x{v)\TL 
'o 



V + UJ 



(3) 



Here A denotes the coupling constant, A^(0) corresponds 
to the density of states at the chemical potential and lOc 
is an upper cutoff energy for the interaction. 

In the limit P ^ where F(jQ_r^o(w) = 5{ijj — ujo) 
Eqs. H2I3|) describe the self-energy as arising from the 
coupling to a single oscillator mode. In this case one 
obtains the well-known step-function behavios^Si of Im 
E(tj) = l/2\^N{Q)Q{\uj\ -ujo). Let us consider for illus- 
trative purposes at first a rectangular distribution func- 
tion i^;^p_r(w) with halfwidth P and height 1/(2P) for 
which the self-energy can be calculated analytically. As a 
result for finite but small P < the step in Im E(cl;) will 
evolve into a linear increase between ljq — T and cjo + P. 
Finally for T — uj^ the MFL limit is reached where Im 
E(a;) ~ w for < w < cjQ + P. 

Furthermore we evaluate the quasiparticle weight Z — 
(1 — {dTi/duj)Lo=Q)~^ which is given by 



Z 



1 H In- 



P 



(4) 



and vanishes logarithmically when P approaches the 
mean mode frequency ojq. This result can be connected 
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FIG. 1: Im E/fcT as a function of u/kT for 48 K (solid line, 
diamonds), 90 K (dashed line, triangles) and 300 K (dotted- 
dashed line). Symbols are experimental data taken from Fig. 
4 of Ref. Ig. Inset: r(T) as obtained from a comparison 
between eq. ||1J and experimental data of Ref. (circles) . 



with recent ARPES experiments where the quasiparticle 
weight Z as a function of temperature has been extracted 
from the line shape at (tt, O).^^- In fact it was found that 
Z drops significantly as temperature approaches Tc with 
a small but nonzero Z also above Tc- From Eq. Q we 
thus conclude that the bosonic spectrum is governed by 
a distribution function F^gx{uj) with a strongly temper- 
ature dependent halfwidth P = P(T). 

We have extracted T{T) by comparing eq. Q with the 
experimentally determined quasiparticle weight (Fig. 4b 
in Ref. IT3) . The result is displayed in the inset of Fig. 1. 
This plot shows that in the SC state P(T) continuously 
increases and reaches the MFL limit (i.e. P « loq) for 

With the extracted tanh— like form for P(T) we now 
evaluate the temperature and frequency dependence of 
Im E for parameters which also fit the real part (see be- 
low). In Fig.l we compare our results with ARPES data 
on an optimally doped Bi2212 sample^ which support 
a scaling of Im E as a function of u/kT characteristic 
for a MFL. However, from the remarkable agreement of 
our fit with experiment as shown in Fig. 1 we conclude 
that the presence of a collective mode in the 48i4r data 
cannot be excluded. We will see below that this is also 
consistent with the observed temperature dependence of 
Re E as required by the Kramers-Kronig transformation. 
It is worth noting that P(T) was obtained from a fit to 
the quasiparticle weight at (tt, 0) whereas the experimen- 
tal data shown in Fig. 1 are extracted from momentum 
distribution curves along the nodal direction. For large 
frequencies Im E \^N{Q) and thus in our scaled repre- 
sentation Fig. 1 the saturation decreases with increasing 
temperature. Note that this saturation is due to the up- 
per cutoff frequency of our bosonic spectrum at frequency 

UJQ + P. 

Fig. 2a displays Re E for parameters appropriate to 
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fit the data extracted from ARPES experiments on an 
optimally doped Bi2212 samplei^ Note that parameters 
are the same for both temperatures except the broaden- 
ing r which in the normal state is five times the value of 
that in the SC state. From our fit we thus again conclude 
that the increase of temperature alone is not sufficient in 
order to explain the change of Re S from being signifi- 
cantly peaked at T = 70K to the much broader structure 
at r = 130K. Instead the transition from T < Tc to 
T > Tj, is accompanied by the transformation of the self- 
energy from mode-like to almost MFL behavior due to 
the broadening of the underlying bosonic spectrum Eq. 



experiments on underdoped Bi2212 materials (cf. Fig. 
2c of Ref. i3) clearly support the existence of a 'break' 
also in the nodal direction corresponding to the situation 
shown in Fig. 2b. Increasing doping or (and) increas- 
ing temperature leads to a decreased change of slope (i.e. 
less pronounced kink) around a characteristic binding en- 
ergy (cf. Figs. lb,e of Ref. 0). Comparing with the 
spectra shown in Fig. 2b-2d all these data suggest that 
the charge carriers are coupled to a bosonic spectrum 
which is rather sharply peaked well below Tc but grad- 
ually broadens with increasing temperature. In addition 
from the doping dependence of the kink we conclude that 
the minimum F is smaller in the underdoped regime. 
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FIG, 2: a) Re ^{u>) compared with experimental data from 
Ref, @ for 70 K (solid line, circles) and 130 K (dashed 
line, squares). Parameters: ujo — 55meV, uic ~ 400meV , 
A^iV(O) = 20meV, r{70K) = 15meV, r(130A') = 75meV. 
(b-d) Contour plot of spectral functions for increasing values 
of P and temperature showing the gradual weakening of the 
kink feature. 

In Figs, 2c-2d we show a contour plot of the spectral 
function Ak{Lu) = l/7rlml/(ti; — Efc — S(w)) calculated for 
the two values of F as determined in Fig, 2a, The un- 
derlying bare dispersion is linearized around the Fermi 
momentum ~ (/c — kp) since experimental data for Re 
S in Fig, 2a have been obtained by substracting the mea- 
sured nodal dispersion from a noninteracting, linearized 
dispersion, ^•^ 

In addition Fig, 2b displays the spectral function for 
a rather narrow boson distribution (F = ImeV) at tem- 
perature T=OK, In this case the contour plot (Fig, 2b) 
reveals the occurence of a break in the dispersion sepa- 
rating a strongly mass enhanced low energy part from a 
broad dispersive high energy feature. With increasing F 
(Figs, 2c, 2d) this break gradually changes into a kink 
which in the MFL limit (F ^ ujq) reduces to a slight 
change of slope in the dispersion. 

As already discussed in the introduction the kink in 
the dispersion along the nodal direction is now well es- 
tablished in a large class of cupratesi Moreover ARPES 



III. BILAYER SPLITTING 

Having analyzed the temperature and frequency be- 
havior of the bosonic spectrum we now proceed by in- 
vestigating its influence on the experimentally observed 
bilayer splitting in the SC and normal state of Bi2212 
materials, respectively. 



Formalism 



Starting point is the following hamiltonian 
H - ^[ek-A^]c,t(k)c,,,(k)-f J2 ix(k)4,(k)c,-.(k) 



+ ^ A(k) [cl (k)4^ (-k) + ca (-k)c,,T (k)] 



(5) 



where — 1,2 label the planes of the bilayer sys- 
tem. We restrict on an in-plane d-wave SC gap A(k) = 
(Ao/2)[cos(fc2:) — cos(fcj^)] and the interlayer hopping 
is parametrized as t±{k) ~ {t^j_/A)[cos{kx) — cos(fcj^)]^ 
in agreement with LDA24 calculations and ARPES 
experiments^^ The bare energy spectrum Ck in each 
plane is a tight-binding dispersion including up to fifth- 
neighbor hopping.-^ A comparison of the corresponding 
Fermi surface (FS) with ARPES data from Ref, jli is 
shown in the inset of Fig, 3c for an interplane matrix 
element t']_ ~ AQmeV . Diagonalization yields the four 
eigenvalues ±Ei^2{k) ^ ±y/{ek - M ± t±{k)y + A(fc)2, 

The system Eq. ^ is now coupled to ICDW type 
fiuctuations via the action 



S 



E 

9 



i,3 



dn I dT2y ,pI{ti) x'H (ti - r2 ) pi , (ra ) 

(6) 

where the susceptibility matrix contains both intra- 
(xl = x*^-') and interlayer (x^ = Xq'^"') contributions. 
Following the approach in Ref. we consider an in- 
plane susceptibility which is factorized into an lo- and 
q-dependent part, i.e. xl(*'^) = •^|W^(*^)>^(q)- Here, A|| 
denotes the intralayer coupling constant and W{iLo) = 



4 



/ dvB{u)2u/{u 



is determined by the bosonic 



spectrum Eq. ^for which we now consider a lorentzian 
shaped distribution function F^^^y{'-^)- Jil) contains the 
charge-charge correlations which are enhanced at the four 
equivalent critical wave vectors (±g^, igp. 



^(q) 



4 



E 



7 



72 + 2 - cos{q^ - q%) ~ cos{qy - q") ' 



(7) 

N is a suitable normalization factor introduced to keep 
the total scattering strenght constant while varying 7. 
Note that for the chosen sign convention in Eq. 10 this 
kind of scattering leads to an effective attraction (in the 
static limit) between holes in the same plane. However, 
for the interlayer susceptibility the situation is different 
since the Coulomb repulsion perpendicular to the lay- 
ers is only weakly screened and it is thus more likely to 
assume an effective repulsive interaction between holes 
in adjacent layers. For simplicity we therefore take in 
the following Xq — ^xf with A|| = Aj^ = A in order 
to describe the physical situation where charge enhanced 
regions in one layer correspond to charge depleted ar- 
eas in the other layer. In this regard the situation is 
analogous to AF correlated spin fluctuations between the 
planes^ when the spin-exchange potential couples exclu- 
sively even and odd electron states. We therefore con- 
sider our approach as complementary to models which 
are based on the magnetic resonance as the relevant 
collective mode. Preliminary investigations on charge- 
instabilities in a bilayer Hubbard-Holstein model^^ also 
give evidence for the presence of ICDW fluctuations in 
the odd channel. 

The leading order one-loop contributions to the intra- 
{a — \\) and interlayer (a =_L) self-energies read as 



q.ip 



q,iuj-ip)T^ (8) 



which in turn allows for the calculation of G via 



G — Co -t- Go ^G 



(9) 



where the unperturbed Matsubara Greens function and 
self-energy are 4x4 matrices and are given by 




Gq Gq" 



Finally by inverting the 4x4 matrix equation (PJ 
the spectral function can be extracted from Ak{u!) = 
/m[Gii(fc,L^) + G33(fc,c^)]. 



B. Results and Discussion 

In Fig. 3 we present results for energy distribution 
curves (EDC) I{k,uj) = f{uj)A{k,Lo) (/(w) being the 



Fermi function) in the normal (T=90K) and SC (T=10K) 
state and wave- vectors are indicated in the inset of Fig. 
3c. Parameters (cf. caption to Fig. 3) were chosen in or- 
der to reproduce the ARPES line shape features in Rcf. 
(Ij. Consider first the normal state spectra shown in Fig. 
3a. Due to the large value of F in the underlying boson 
spectrum the bonding band (BB) gives rise to a broad 
high energy feature in the EDC at k = (tt, 0). Together 
with the antibonding band (AB) quasiparticle peak the 
EDC thus displays a normal state peak-dip hump struc- 
ture. The coupling to the fiat boson spectrum naturally 
also affects the peak-hump separation at k = (tt, 0) which 
therefore is slightly larger (« QOmeV) than the bare bi- 
layer splitting (2^^ = 80meT^). The MFL structure of 
both AB and BB spectral line shape is resolved in Fig. 
3c. Upon scanning from (tt, 0) towards (tt, tt) first the 
AB peak goes through the chemical potential while the 
BB peak shifts to lower energy and sharpens due to the 
Im S ~ a; behavior. Finally between nS and n9 also the 
BB peak crosses the Fermi energy. The AB and BB peak 
positions of the EDC curves from Fig. 3a are compared 
in Fig. 4 with ARPES data from Ref. d 
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FIG. 3: ARPES spectra calculated for (a) normal state 
(T=90 K, Ao^OmeV, r=100meV), (b) SC state (T=10 K, 
Ao=20meV, r=10meV). Spectra are taken along (7r,0) 
(7r,0.247r), and labeled from to 9 as shown in the inset of 
(c). The inset displays the calculated FS in comparison with 
data from Ref. (c) shows selected spectra from (a) and (b) , 
where dashed and dotted curves indicate BB and AB states 
respectively. Parameters: A=0.095eV, qc = 7r/4, a;o=40meV. 



In the SC state the narrow bosonic distribution (F — 
IQmeV) induces a strong mass enhancement for quasi- 
particles in the binding energy range up to the mean 
mode frequency ojq. Both AB and BB line shapes display 
the standard peak-dip hump structure where the hump 
is more pronounced for the BB (cf. curve s4 in Fig. 3c). 
The stronger renormalization of the BB band is due to 
the fact that (for our choice Xq — ~Xq) self energies 
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of AB and BB bands 

(10) 

are determined by fermionic excitations of BB and AB 
bands, respectively. Since the van-Hove singularity of the 
BB band is well separated from the chemical potential 
this results in a small AB self-energy and thus a rela- 
tively weak AB hump feature. In contrast the closeness 
of the AB van-Hove singularity to fi induces large BB 
self-energy effects and therefore a more pronounced BB 
hump structure. 




120 I ' ' ' ' ' 1 
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FIG. 4: Dispersion extracted from Fig. 3 (open symbols); 
Filled symbols indicate experimental data from Ref. Q. Cir- 
cles and triangles correspond to BB and AB normal state 
peaks, respectively; Diamonds and squares refer to BB and 
AB quasiparticle peaks in the SC state. 

The SC low energy 'quasiparticle' peaks are confined 
to energies smaller than loq and the corresponding quasi- 
particle weight decreases with increasing binding energy. 
As a consequence the quasiparticle peak of the bond- 
ing band is strongly suppressed at k = (tt, 0) (curve sO 
in Fig. 3b) and reemerges when the bare energy of the 
bonding band approaches luq (curve s2 in Fig. 3b). In 
Fig. 4 where the low energy peak positions are shown as 
a function of kx the dispersion of the BB therefore seems 
to be interrupted around (tt, 0) in agreement with the ex- 
perimental data of Ref. Q. Moreover, Fig. 4 displays one 
of the main results of this paper. Whereas in the nor- 
mal state the bilayer splitting is around 'dQmeV the low 
energy peaks in the SC state are separated by « 20meV 
only. As already discussed above the decreased bilayer 
splitting in the SC state is only due to the sharpening of 
the collective mode at low temperatures and not influ- 



enced by any other parameter in our calculation. 
IV. CONCLUSION 

To conclude, we have shown that the self-energy as 
extracted from ARPES experiments can be consistently 
described within a model where the charge carriers are 
coupled to a bosonic distribution centered around some 
mean frequency uq. The width F of this spectrum turns 
out to be strongly temperature dependent as can be 
concluded from the temperature dependent quasiparti- 
cle weight (cf. inset to Fig. 1), the evolution of the 
kink along the node (Figs. 2b-2d) and the correspond- 
ing change of Re S in Fig. 2a with temperature. Within 
the ICDW scenario the associated scattering can account 
for the experimentally observed^ appearance of an addi- 
tional peak-dip hump structure below Tc near the (tt, 0) 
points of the BZ (curve s4 in Fig. 3c). Moreover, the 
strong QP mass enhancement below Tc due to the small 
value of F leads to an apparent reduction of the bilayer 
splitting and a suppression of BB spectral weight around 
the (7r,0) points (Fig. 4). 

Note that phonons cannot account for the tempera- 
ture dependence of the collective mode although their 
linewidth may undergoe some sharpening upon enter- 
ing the SC state. For temperatures well below the 
phonon frequency one can estimate the linewidth^S from 
F « ^A^(0)AijJo which yields F << ojq also for large cou- 
pling A « 1...2. This is not enough to account for the 
strong damping as discussed above. 

On the other hand the temperature dependent width 
of the bosonic spectrum can at least be made plausi- 
ble within the ICDW scenario. At low temperatures the 
large intensity CDW mode at energy luq and wave vecors 
(±gc, ifZc) is part of a band of excitations which for q — > 
corresponds to the phason mode of the system. It is 
knowr^ that the effective mass of the phason strongly 
increases with temperature due to the decreasing CDW 
gap inducing a significant broadening of the bandwidth of 
charge excitations. In fact, preliminary results obtained 
for the stripe phase of the Hubbard model^^ near the 
CDW transition have already revealed a rather continu- 
ous frequency spectrum in the charge susceptibility cen- 
tered around Qc corresponding to the stripe periodicity. 
However, a more detailed investigation of this tempera- 
ture dependence remains an interesting issue for future 
work. 
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